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Abstract. We prove the existence of a wide class of solutions to the isentropic relativistic Euler equa- 
tions in 2 spacetime dimensions with an equation of state of the form p = Kp^ that have a fluid vacuum 
boundary. Near the fluid vacuum boundary, the sound speeds for these solutions are monotonically 
^ ' decreasing, approaching zero where the density vanishes. Moreover, the fluid acceleration is finite and 

^ , bounded away from zero as the fluid vacuum boundary is approached. The existence results of this 

article also generalize in a straightforward manner to equations of state of the form p = Kp i with 
! 7 > 0. 

(N 

O^l 1. Introduction 

U 



(N 



Letting 

g = g^^dx^dx^ /i, = 0, 1 
denote a flat Lorentz metric of signature (+,— ) on a 2 dimensional manifold M C M^, the relativistic 
^ \ Euler equations are 

\0 '. V^.T^" = (1.1) 

, where 

T^'•' = {p + pjv^'v" - pg"" and gf.^v^v'' = 1. 

Here, is the fluid two- velocity, p the proper energy density of the fluid, and p the pressure. Projecting 
I (|l.ip into the subspaces parallel and orthogonal to yields the following well known form of the Euler 

equations 

v^V^p + (p + p)V^w'' = 0, (1.2) 
{p + p)v^V^v'' -h'"'y^p = (1.3) 
' where 

C3 . ^ g^^ - Vfj,v^ {Vf_, = gfj.uv'^) (1.4) 

is the induced negative definite metric on the subspace orthogonal to v^. In this article, we will be 
primarily concerned with fiuids with an insentropic equation of state of the form 

p^Kp"^ K,j>0, (1.5) 

and in particular 7=1. 

The main aim of this article is to construct solutions to the Euler equations that contain a vacuum 
region; that is solutions for which p vanishes, and the fluid vacuum boundary is located where p — 0. 
The problem of existence of solutions to the Euler equations with a vacuum region for equations of state 
for which the pressure and proper energy density simultaneously vanish has been studied by a number of 
authors. The main difficulty in establishing existence is due to the fact that the Euler equations p.2p - (|1.3|) 
become degenerate as the fiuid vacuum boundary is approached. This is due to the vanishing of the the 
quantity p + p at the fiuid vacuum boundary. Because of this, standard techniques from hyperbolic PDE 
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theory do not apply. The first general existence result where this problem was overcome is in |16j . In this 
work, the existence of gravitating, non-relativistic fluid bodies with compactly supported densities in 4 
spacetime dimensions was established (see P^j and ^15 ^ for the relativistic case); the technique employed 
involved a special choice of variables to regularize the Euler equations, and works in all dimensions with 
or without coupling to gravity. However, the type of fluid solutions obtained by this method have, when 
coupled to gravity, freely falling boundaries, and hence, do not include static or expanding fluid bodies. 
This is because the fluid acceleration 

a^'^v^V^.v'' (1.6) 

of the solutions from jl6) vanish at the fluid vacuum boundary, and so there is no outward force to 
counteract the gravitation force which leads to collapse of the body. Consequently, to have a class of 
solutions which include the full physical range including static and expanding fluid bodies, it is necessary 
to establish the existence of solutions to the Euler equations for which the acceleration (|1.6p is non-zero 
at the fluid vacuum boundary. Due to the finite propagation speed of the Euler equations, it is enough 
to consider the existence problem in a neighborhood of the fiuid vacuum boundary. Away from the 
boundary where p + p > 0, standard symmetric hyperbolic techniques can be used. 

In this article, we establish the existence of a wide class of solutions in 2 spacetime dimensions for 
which the acceleration is non-zero at the boundary. The approach we take to establish existence relies 
on geometric arguments and is inspired by the Lagrangian formulation of the relativistic Euler equations 
in [8], although we use a different formulation of the relativistic Euler equations due to Frauendiener 
and Walton [71 HT] . In the Frauendiener- Walton formulation, the normalized two- velocity and proper 
energy density p are combined into a single vector . Due to the vector nature of this formulation, 
it is possible to exploit the full diffeormophism freedom available to fix the coordinates, and reduce the 
problem to that of two linear wave equations one of which is regular while the other singular. Existence 
for the regular wave equation is handled by standard PDE techniques while existence for the singular wave 
equation is obtained by standard results on abstract wave equations. Although the analysis in this article 
applies only to 2 spacetime dimensions, some of the arguments presented here do have a counterpart in 
higher dimensions. However, we will not discuss these results here. 

In 2 spacetime dimensions, the existence of solutions to the non-relativistic Euler equations with non- 
zero acceleration at the fiuid vacuum boundary has been previously established in the remarkable articles 
[5l[Tl]. In these articles, existence is proved using non-standard energy estimates combined with suitable 
approximation techniques. The arguments used are technical, involved, highly original, and quite different 
from one another. We also note that the results in [Sj [11] have been recently extended to 4 spacetime 
dimensions [51 [T2] (see also [1]). At the moment, the relationship between the solutions of this article 
and those of [3 [11] is not clear because the solutions presented here are relativistic while those from 
[5] 111) are non-relativistic. In order to understand the relationship, it would be necessary to take the 
non-relativistic (i.e. c — > oo) limit of our solutions. We will not do this here, but leave the analysis of 
this limit to a separate article. 

Overview of the article: In Section[2l we review the Frauendiener- Walton formulation of the isentropic 
Euler equations, and we use this formulation to show, in Section [3l that the Euler equations are equivalent 
to the existence of a suitably normalized, commuting set of vector fields. This commuting set of vector 
fields allows us to introduce adapted coordinates for the Lorentzian metric. As detailed in Section 14.11 
the fiatness of the Lorentzian metric then implies the equivalence of the Euler equations with a non-linear 
scalar wave equation. Conformal coordinates are introduced in Section 14.21 which reduce, as described in 
Section l473l the non-linear wave equation to a linear one. In Section |4^ it is shown that the existence of 
conformal coordinates is equivalent to the existence of suitable solutions for a linear wave equation with 
singular coefficients. An existence and regularity theory for this type of wave equation is developed in the 
Appendix [Aj The main result of this paper is contained in Section 14.51 where the existence of solutions 
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to the Euler equations with non-zero acceleration at the vacuum boundary is established, see Theorem 
14.81 for a precise statement. Finally, in Section [SJ a class of exact solutions with non-zero acceleration at 
the vacuum boundary are described. 

2. The Frauendiener- Walton formulation of the Euler equations 

In [71 [21] , Frauendiener and Walton independently showed that the isentropic Euler equations for a 
perfect fluid with an equation of state of the form p — p{p) can be written asQ 

A^.-'V^w" = 0, (2.1) 

where w'^ is a timelike vector field with norm 

= w^w'' > {Wf^, = gi_,uw''), 

and 

A^,^ = + ^) ^^^^ - ^>t^ - - (2.2) 

We will refer to these equations as the Euler- Frauendiener- Walton (EFW) equations. 
In the Frauendiener- Walton formulation, is a function of 

w 

where 

w ~ \fuP' . 

An explicit formula for can be calculated in the following fashion (see [7] for more details). First, the 
pressure p — p{(^) is determined implicity by the equation 

C = CoT(MC)), (2.3) 

where 

T(p) = exp ( r --^) (2.4) 



(2.5) 



pip) + p 

is the Lichnerowicz index of the fluid. From this, can be calculated using the formula 

1 ^ fa'io 

\ /(C) 

where 

/(C) = CV(C)- 

Additionally, the proper energy density p and two-velocity can be recovered from 

P = PiC)-^ and w^^C^^". (2-6) 
C 

As shown in [7], and also , the triple u'*} determined from (|2.ip . (|2.3p . and (|2.6I) satisfy the 

relativistic Euler equations (|1.2l) - (|1.3p . 

For the equations of states (|1.5|) . it is possible to explicitly determine the functional form of ~ s'^iC)- 
To see this, we observe that the Lichnerowicz index (|2.4p is given by 



^It is important to note that Frauendiener and Walton use opposite signature conventions for the metric g and different 
notation for the fluid variables. In this article, we use the notation and signature conventions of Frauendiener [J- 
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From this expression, it is clear 

( 1 \ 7+1 

solves (|2.3p . and hence determines the pressure as a lunction of (. Without loss of generality, we set 
Co = 1 which gives 



P(C) = ;^(C^-1) 
Substituting this expression into (j2.5p then yields 



.^ = ^(C^-1). (2.7) 

7 

For latter use, we note that by using (|1.3p the fluid acceleration can be written as 

a" = — ■ — V^p. 

P + P 

Since this vector is orthogonal to v'^, i.e. v^a'^ — 0, we can use the metric /i^^ to calculate the length of 
a'^ via the formula 



\a\h V^V^^ = \/ -7r7-:7;2 ^mPV^p . (2.8) 



Recalling that the square of the sound speed is given by 



2_dp _K{-r + 1) 



dp 7 



we can write (12.81) as 



^ -V^s^V,s^ = —-^\ds^\h. (2.10) 



This shows that the acceleration will be non-zero at the fluid vacuum boundary if and only if | ds^|/i > 
there, since p vanishes at the boundary. It is also clear from this formula that the fluid acceleration can 
only be non zero at the fluid vacuum boundary if the gradient of s"^ does not vanish there. Since the 
proper energy density and hence the sound speed is decreasing to zero as the fluid vacuum boundary is 
approached, it follows that for such solutions, must be monotonically decreasing in the neighborhood 
of the boundary. 

3. A FRAME FORMULATION FOR THE EFW EQUATIONS 

Since the metric g — g^^^dx^dx'^ is flat, we begin by introducing global Minkowskian coordinatef0 
[x^) = {x^,x^) for which 

{9,u)^[l (3.1) 



"^We note that this step of introducing a global Minkowskian set of coordinates is not necessary. We include it in order 
to make the following arguments more accessible to readers who may not be familiar with a more abstract approach to 
fixing coordinates. A more abstract approach becomes particularly important when dealing with non-flat metric that arise 
when coupling the fluid to gravity. 
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With this choice of coordinates, the Euler equations (|1.2p - ()1.3p can be written exphcitly in terms of the 
proper energy density and a function -y(a;°,x^) that satisfy 



with the two-velocity given by 



K) 



-A pi 

1 

1 



1 — 



r^oP + 1 + 



dip]=Q 



Using (|2.7p and (|2.9p . a short calculation shows that 

1 



w 



from which we obtain 



by (|3^ and ([2S1)- Letting 



(1 + 
1 



:(l,w(a;0,a;i)) 



and 









2 








/ ^2 








9 








/ 



1 \ 



0\2 



W 



3+ ^ 



1 \ {w^y 



1 \ (u;i) 



1\2 



3 + - 



1 \ {w^y 



where 

and (see (^71) 1 



U;^ = («,0)2_(^l)2 



7 + 1 



7 \(w2) 2(7 + 1) 

we can also write the the EFW equations (I2.ip explicitly as 



(3.2) 



(3.3) 



In order to derive a formulation of the Euler equations that is suitable to analyze the limit s'^ \, 0, we 
use a particular frame formulation of the EFW equations. First, we set 

where 

d 



(3.4) 
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and choose 

ei = e^d^ 

orthogonal to cq so that the frame metric 

9ij = g{et,ej) = Spi^efeJ' 

satisfies 

501 = 501 = 0, (3.5) 

and 

500 = 5(eo,eo) = w^- (3.6) 
At the moment, we leave the length of ei unspecified. The freedom to fix the length of ei will be used 
below. From, p.ip and p.3p . it is clear that 

{et)=q{x'',x'){v{x",x'),l) (3.7) 

where qix'^, a:^) is a function to be determined. 
Next, we denote the coframe by 

= dldx'^ (0^6^ = <5}), 
and let uji'^j denote the connection coefficients so that 

We also define the connection 1-forms ui'^j in the standard fashion 

and set 

^kj = gu^^ j — ^ikjO^, 

where 

t^tkj = gu^i] = i^ik] = 5(Veiej, Efc). (3.8) 

Letting 

A*^' = 5 Ve^e^^^^^.^ 
a short calculation using p.2p . p.4p . p.Sp . and p.6p shows that 

A*^"^- = ^3 + 1^ - Slg^" - g'Hl - g'^^. (3.9) 

V s / 500 

Moreover, it follows from p.4p and p.Sp that 

Wfcjo = e^V^w^eJ^i,^. (3.10) 

Together, (13. 9p . p.lOp and the invertibility of g*^ and show that the EFW equations (|2.ip are 
equivalent to 

A'^'ujkjo = 0. (3.11) 
Using (1^ . ((Xa . dMI, and (jX^ . the i = and i = 1 components of (|XTT|) are 

1^000 1^110 = and cjioo + ujoio = 0, (3-12) 

511 



respectively, where 

2 7+1 



7 V5oo"(^+i' 



1 . (3.13) 



Also, due to the connection w'j being metric, the frame metric satisfies (see [2j Ch.V,§B]) 

dgjk — ujjk + ^kj 
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or equivalently 

eiigjk) = Wjjfc + i^ikj, (3-14) 
and it follows immediately from p.Sp that 

Wool + t^oio 0. (3.15) 

Since the connection aj'^j is torsion freeH, it satisfies the following Cartan structure equation (see [U 
Ch.V,§B]) 

d0' + oj'j A0^ =0, 

or equivalently 

[eo,ei] = (ujo'^i - i^i''o)ek- 

We rewrite this as follows: 

[eo, ei] —{ujQio + wioo)eo + — (won - wiio)ei by (|33|) and p.lSp 

.900 .911 

= — ( Won l^wooo ) ei by p.l2p 

511 V s-ffoo / 

^ ' eo(ffii) - #^eo(5oo) ) ei. by (EH. (3.16) 



2011 V s^ffoo 

Defining a function by 



m = — ^ 



^2FfTT) ( \ - 1 

a short calculation shows that F(f ) satisfies 

= pi 



7 

and hence 

i^'(3oo 

by p.l3p . But this implies that 



j^(.9oo) 
2500S 



2 



eo(ln(F(goo)^) = ^^60(500), 



s^5oo 



which we can in turn use to write (13. 16^ as 



[eo,ei] = ieo(^ln(^^^jjei. (3.17) 

We now use the freedom to fix the length of ei by setting 

511 = -Figoof. (3.18) 

By p.l7p . we arrive at the equivalence of the EFW equations with the vanishing of the following Lie 
bracket 

[eo,ei]=0. (3.19) 

We also note that if we set 



■^This follows by virtue of the connection u'^j being metric. 
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then ()3.20p can be solved for goo to give 

w'=goo = {l-uf^-'+^\ (3.21) 
which allows us, using (|2.7p . to write the square of the sound speed as 

(3.22) 



2 _ 7+ 1 u 



7 1 — u 

Using this and the normalization p.lSp condition, we see from p.7p that 

1 1 



7±1 +^2 
7 / 



and hence, that 

i4:T/9(a;0,xi)^l - i;(2;0,xi)2 ' 

by (USD. 

To summarize, the main result of this section is that the EFW equations are equivalent to the vanishing 
of the Lie bracket p.l9p which in components reads 

where the frame components are given by 

(eo) = , V, , (l.^;r.T°.xM). (3.23) 

4. Existence of solutions to the Euler equations 

4.1. A wave equation formulation of the EFW equations. The vanishing of the Lie bracket p.l9p 
implies the existence of coordinates which, at least locally, trivialize the vector fields cq and ei. To 
construct these coordinates, we let 

j-,,,(x",ii)- (j:o,(iO,xi),j-i,(xO,ii)) i = o,i 

denote the flow maps of the vector fields {i — 0, 1), that is J-i^r is the unique solution to the initial 
value problem 

^Kr^i\ i') = er(j-.,.(x°,ii)) M = 0, 1, 
J-,,o(£°,£^) = (i°,xi), 

where the ef are given explicitly by the formulas p.23p - p.24p . By translating, if necessary, we can 
assume the origin [x^^x"^) = (0,0) is in the domain on which the vector fields ef are defined. We then 
introduce a change of coordinates by the formula 

{x\x^) = -^{x\x^)^F^ /^ioO J-i,£l_e(0,0). (4.1) 



where c is a constant^. 

Since the vector fields w — eo, and ei commute, if we define 

w — Co — 5'*eo and ei = ^*ei 



*The constant c is chosen so that the point hm^i^y ^i—cCi 0) lies on the vacuum boundary where p = 0. 
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then 



where 



w = eo — \l ^- — do and ei = di, 



d 



(4.2) 



As we show below (see Remark 14. 7p , the vacuum boundary where p vanishes is contained in the set 
= 0. Since the two-velocity is given by v'' = (w)~^w^, equation (|4.2p shows that the vacuum 

boundary moves with the fluid as expected. 
Next, defining 

g — ^* g = (jp^ijdx^dx^ and ii = 
it follows directly from g^), (jST^ . (|3:20)) . and (|33T|) thalQ 

-yH - 7-}'l2(7 + l) 1 



5 



7+1 



-dx dx di 



Computing the Ricci scalar of this metric, we find that 



R = 



G(w)i/2 



2(7+1)4 



1 



(4.3) 



(4.4) 



where 



G(m) = (1 -U)2(T+1). 

In order to simplify the following calculation, we will from this point on assume that 7 = 1, and define 



z — arcsin(— 1 + 2u) + 



(4.5) 



which can be inverted to give 

u = i(l - cos(i)). (4.6) 

Remark 4.1. To extend the analysis to 7 > 0, the appropriate z variable that replaces (|4.5p can be 
obtained by solving the initial value problem 

1/2 

: z(0) = 



dz 



G'{u) 



du \ 2(7+ l){tG(w), 
for u > 0. A solution to this initial value problem yields the identity 



1 



dz 



iPG'{u) dz 



u-r+^Giuy/^ du V 2(7+ l)G(?i)i/2 
which allows a similar analysis as in the 7 = 1 case to be used. 

In terms of the z variable, the metric (|4.3p and the Ricci scalar (|4.4p become 

4 



and 



R = 



^ (l + cos(i))^ ^.o^.o 
^ 32 

(1 — cos(i;)) 



rdx^dx^ , 



2(1 +cos(z))2 



do 



sin(z) 



(1 - cos(z))2 
doz]-dJ''''^'^' 



diz 



(4.7) 



(4.8) 



^The metric g = 'l'*g is just the original Minkowski metric g (see {XT}) expressed in the (£'') coordinates defined by 
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respectively. But the metric (|4.7p is and so (|4.8|) and the above arguments show that the EFW 

equations (|2.2|) are equivalent to the wave equation 

do ( -4^do^) - 5i ( '^d^z] = 0. (4.9) 



ySin(z)-^ y \ " / 

Remark 4.2. It is well known that in 1 + 1 dimensiontQ the Euler equations can be reduced to a quasi-linear 
scalar wave equations of the form (see 3 or [19^ for details) 

V^(ff(|V^|2) V'V) = (|V^|2 = g'^'-d.ipd.ip), 

where H{-) is a particular function determined by the equation of state, and the proper energy density p 
and two-velocity can be recovered via the formulas 

In light of this, it is perhaps not surprising that we are able to reduce the Euler equations to a quasi-linear 
scalar wave equation. What is new here is that proper energy density is a function of the scalar z, and it 
vanishes where z vanishes. As will be shown below, this makes the wave equation ()4.9p particularly well 
suited for analyzing the vacuum boundary problem. 

Defining an auxiliary Lorentzian metric A bjQ 

A = X.^dx^dx" = — ^-^dx°dx° —dx^dx^, (4.10) 

8 sin(z)'' 

equation ()4.9p is easily seen to be equivalent to 

a-^z = -i=4 f v^A^'^a.A = 0, (4.11) 



'-A 

where 

V^= y/-det(V) and X^"" = 
4.2. Conformal coordinates. We introduce a second change of coordinates 

{x°,x^) = ^{x\ x^) = (^(S°, x^), 4>{x°,x^)) (4.12) 

and let 

z = (4.13) 
Pulling back the metric ()4.10p using the diffeomorphism ()4.12p . we find, using ()4.13|) . that 

A = 4>*A = A^^dx^dr, (4.14) 

where 



(V) 



^ /x(9oV')^ - -(9o<^)^ p,doil}diii - -do(t)di(j) ' 

^ _ _ f _ _ 

>/i9oV'i9iV' dQ(j)di(f> nidiip)'^ {di(j)) 

_ _ sin^(z) 



''We recall that the Ricci scalar is a coordinate invariant which can be stated as R = R o ^ . Since g is the Minkowski 
metric and consequently flat, its Ricci scalar vanishes. By the coordinate invariance, the Ricci scalar of g must also vanish. 

'7 

In fact, this is true for irrotational fluids in any dimension. 
^The metric A^,y is conformal to the acoustical metric dfu, = s^VuVv + hf^i, with the relation between the two given by 
_ siii^(2) \ 
"M" - 2(l-cos(z))-^/"^- 
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and 

. _d_ 

We fix the choice of coordinates by requiring that the metric (|4.14|) is conformally flat in these coor- 
dinates. We accomplish this by demanding that and ip satisfy 

fldo-ip = di<f) and = fidi'tjj. (4-15) 

Witli this choice, the metric (|4.14l) can be written as 

A = Hm, (4.16) 

where 

m = dx^dx^ - dx^dx°, (4.17) 

and 

ft 

Also, from the Jacobian of the transformation (|4.12p 

■'♦ = (11 Ij)- 

we see that 

detJ^^n. (4.19) 

4.3. The conformal wave equation. By the conformal invariance of the wave equation in 2 dimen- 
sion^, and the invertibility of the coordinate transformation ()4.12|) on regions where det does not 
vanish, it follows from (|4.16p and (|4.19p that the wave equation (14. lip is equivalent to 

□,ftZ = B^z - dfz = 0. (4.20) 

We will solve this system on the spacetime region 

U = {ix°,x^)\x°,x^ >0} 

with the boundary condition 

z\r=0, 

where 

r^{{x°,o)\x° >o}. 

Since the vanishing of z implies the vanishing of the proper energy density p (see (j2.9p . p.22p . (|4.5p . and 
(|4.6p ). the boundary condition will ensure the the vacuum boundary lies in T. 
We let 

E = {(0,xi)|xi >0} 
denote the initial hypersurface, and and we define the following related sets 

Ut,s ^ {ix°,x^)\0 < x° <T and 0<x^<d}, 
rT = {{x°,0)\0<x° <T}, 

and 

= {(0,x^) |0 < x^ <S}. 



^The conformal invariance follows directly from the identity Oq^z = ^O^z. See Appendix D of |20| for a derivation 
of this identity. 
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= 0, 




= 0, 




= Zq 




= Zl 



Proposition 4.3. Suppose k € N, s > 1/2 + k, zq e H^iT,), zi e H''+'^{T,), and the 

ze^dliZi-2 i = 2,3,...,s 

satisfy the compatibility conditions 

Zi\x^=o^O ^ = 0,l,...,s 

where E and T meet. Then 

(i) there exists a unique solution z £ C^{U) to the initial boundary value problem 

(4.21) 
(4.22) 
(4.23) 
(4.24) 

(ii) for any T,5 > 0, 

z£C\UT,s)nW''^°^{UT,5), 

and 

(iii) if < 2k < diz\^^ < /? and \doz\-Ss < k/4 for some k, 5, /? > , then there exists a Tq such that 

{id,zr-idozr)\^^^>Ky4>o^ 

0<K<9izL < 2/3, 



KX^ < z{x°,x^) < 2Px^ 

for all (5°, 5^) G Uto,5- 

Proof. Statements (i) and (ii) follows from standard linear hyperbolic theory for initial boundary value 
problems. For example, see Theorems 3, 4, 5 and 6 in Section 7.2 of [1]. The first two bounds from 
statement (iii) follow from the bounds on the initial data and the continuous differentiability of z{x'^, x^). 
The last bound in (iii) follows from integrating the second with respect to x^ while using the fact that 
z{x°,0)=0. □ 

Remark 4.4. 

(i) Clearly, we can explicitly solve the initial boundary value problem (|4.2ip - ()4.22p using the well 
known formula 

_ i (zo(x'' + - -zoix" - x') + z,{0 d^) x°<x\0<x' 

'"^ ^ {zo{x° + i^) + -zo{x' - x^) + J^^:^^o m) di) x' <x^,0<x' ' ^^'^^^ 



The point of Proposition 14.31 is that it extends in the obvious manner if we replace U by an 
open set that is bounded by a timelike hypersurface T and a spacelike hypersurface S. In this 
situation, there is no equivalent to the simple formula (j4.25p for solutions to the initial boundary 
value problem (|4.2ip - (|4.22p . As we shall see below, the only property of U that we need is 

{^x = iCx°,^x^) \ X eU and < ^ < 1 } C C/. 

In fact, this can be weakened to the existence of a e > small enough so that 

{^x={^x°,Cx^)\xeUnB^iO) and < ^ < e } C C/, 



where B,{0) = {xeR'^ \ \x\ = V(^°F+(^ < e }■ 
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(ii) It is not difficult to show that the conditions imposed on the initial data in Proposition 14.31 are 
satisfied for a wide class of initial data. For example, initial data of the form 

zo = cy+p{y), zi^q(y) (c = const >0), 

where 
(ii.a) 

Piy) = o(?/*) and q{y) = o{y''^^) as y \ 
if s is an even integer, and 

(ii.b) 

piy) = oiv" ^) and q{y) = o{y' ^) as ?/ \ 
if s is an odd integer, 
satisfy all the conditions on the initial data in Proposition 14.31 

4.4. Fixing the conformal coordinates. Letting denote the Hodge dual operatoiF^ of the metric 
(|4.17l) . we can write (|4.20p as * d(* dz) ~ which implies that the one form * dz is closed. Consequently, 
we get by the proof of the Poincare lemma (see Ch. V §4, Theorem 4.1 of [13]) that 

z{x)=[ Bizi^x) x" + dozi^x) x'^ + zo (zo e M) (4.26) 



satisfies 

dz = *dz, (4.27) 

or equivalently 

Bqz = diz and diz — Bqz. (4.28) 

Evaluating z ai — 0, we see that 

z(0,x^)=[ zi{0,^x^)x^ + Zq. 



Since zi = o(a;^) as X"'^ \ (see Remark l4.4l fii)). it is clear that by choosing zq to be 

zq = — min < 0, inf / zi{0,£^x^) x^ d^> 
I o<si«5Jo J 

we can ensure that 

0<z(0,i^)<(5 forO<x^<(S. 
It then follows from the integrating the first equation of (|4.28p and the bounds on diz from Proposition 
[331 (iii) that 

KX° < z{x°,x^) <2l3x° + S (4.29) 

for all {x^,x^) g UtoM- 

Since dz A dz ((9iz)^ — (^o^)^) dx° A dx^ , the non-vanishing of ((9iz)^ — (9oz)^) on Uto,s (see 
Proposition 14. 3p shows that {dz,dz} forms a basis for the space of one forms at every point of Uto,s- 
This allows us to look for solutions to (|4.15l) that are of the form 

<f> = <j>{z,z) and -ip — iIj{z, z). (4.30) 

To see this, we get from (li??7)l and ([T^ that 

*m d(f) = dz(j)dz + dz(t)dz, (4-31) 

and 

d^ = dz4> dz + dzij} dz. (4.32) 



^'^We recall that the Hodge dual operator on one forms is defined by *m dx^ = y/\m\m>^'^ e^a da;"" where |m| = — det{m^^) 
and e^o- is the completely antisymmetric symbol. In particular, this implies that *m ASp = Ax^ and *m Ax^ = dx'^. 
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Writing (|4.15|) as 
where 



(4.33) 



we see from (|4.3ip . (|4.32p . and (|4.33p that ^ and satisfy 

IJ.{z)dz'ip = dzcf) and fj,{z)dz7p = dz4>- 
From these equations, we then get the wave equation 



5,^</>-M(z)9.-( = 



(4.34) 
(4.35) 



for 



Proposition 4.5. Suppose the initial data {4>\z=aidz4'\z=o) for the wave equation (|4.35p satisfieJP\ 



z=0 



z=0^ 



520(O,Z) > c> 0, 



(4.36) 
(4.37) 



and 



dz(t)iO,z)>0 (4.38) 
/or all < z < tt/2. Then there exists a unique solution (j){z,z) to (|4.35p that can be written as 

^(z,z) = ^(z,z4), (4.39) 

where (j)(t, ^) satisfies 



and 



4> e ([0, oo), C'-'/^O, {^/2f)) n fl ([0, oo), C5-^-(0, (7r/2)4)) 
|<9*^(t,OI + l'^(i,OI<e 



Moreover, there exists a Tq such that 
for all (t,0 e [0,To) X (0,7r/2). 

Proof. This is just a restatement of Theorem IA.4I from the Section lA. 31 of the Appendix. 



(4.40) 

(4.41) 
(4.42) 

□ 



Integrating the first equation in (|4.34p with respect to z, and using the second equation of ()4.34p to 
fix the undetermined function of integration, we find the fohowing expression 



ipiz,z) 



m(^) 



dz(l>{T, z)dT + / -—dz(f>iOX)dC 



M(C) 



for ip. Using (|4.40p . we can write tp as 



32z 



3 /-z 



sin (z) Jo 



sin^(C) 



a,0(O,C^)dC- 



(4.43) 



(4.44) 



^^The spaces tiJ^ are defined in Section lA. 21 of ttie Appendix. 
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Lemma 4.6. Suppose k ~ 3, z is the solution to the wave equation (|4.20p from Proposition \4:.3\ z is a 
defined by (j4.26l) , (f) and cj) are the maps from Proposition 14.51 and ip is given by (j4.44p . Then for Tq and 
Sq small enough, the change of coordinates map ()4.12p 

^ : C/to> : (x) ^ {^{x),^ix)) - z(x)) , </)(z(S), z(x))) 

is well defined and of class on Utq.s, md the Jacobian matrix is given by the formula 

( d,M<x),z{x)) d,^{z{x),z{x)) \ ^ , . 

Mx)^\ ^(z(x)) ^i{z{x)) \[i-T-\ T-T-{ 

\d,(p{z{x),z{x)) d,4>{z{x),m) I 

and satisfies 

detJ^\uTo,5 > 0. 

Furthermore, 

f.x° diZ(ix° ,0}di+zo 



^{x°,0)= (32 J d^^{T,0)dT,0 



Proof. By Proposition 14. 31 fiii) and (|4.29p . we have that 

{z{x), z{x)) e [0, 2(3Ta + S) x (0, 2(36) 

for ah X € Uto,s- Choosing 6 and Tq small enough, it is clear that we can arrange that 

{z{z),z{x)) e [0, To) X (0, 2/3(5) C [0, to) x (0, 7r/2) (4.45) 

for all X G Uto,s- As a consequence, the change of coordinates map (see (|4.12p ) 

(x) = ^(x) = {i'{x),^{x)) (V^(z(x),z(x)),0(z(x),z(x))) (4.46) 

is well defined for all x e UTf,,s 

A short calculation using (|4.18p and (|4.46p shows that 

dzip{z,z) dztp{z,z)\ (Bqz Biz 



'^^ \dz4){z,z) di(j){z,z) ) \doz diz) ' 
Using (|4.28p and (|4.34p , we can write this as 

'ds(j){z,z) dz(p{z,z) 



Taking the determinant gives 

- "<« {iWJ - (^)l «^'-''' - 

It follows from (14.391) and (|4.42p that there exists a positive constant C such that 

dsc^jz, z) 4z3 -4n ^ ^ ^ n 

7-— = —T-rd^n^, z )>C >0 

fi{z) fi{z) 

for all {z,z) e [0,To) x (0,7r/2). This, in turn, implies via (|4.45l) that 

dz(t){z{x),z{x)) 4:Z^{x) 



H{z{x)) m(^W) 

for all X e Uto,5- 



d^(l){z{x),z'^{x)) > C > (4.48) 



16 



T.A. OLIYNYK 



Next, we observe that 



d^(j){z,z) 



dt4>{z, z^) 



< 



fj,{z) 



for all (z, z) e [0, To) x (0, 7r/2) by (|4:39| and (lOTj) . This mequality, with the help of (1445)) and the fact 
that linij^o = 0, shows that by choosing S small enough, we can arrange that 



d:,4>{z{x),z{x)) 



< 



C_ 

71 



(4.49) 



[i{z(x)) 
for all X G Uto,5- 

The two inequalities (|4.48l) and (|4.49l) together with the bound on (diz)'^ — (doz)"^ from Proposition 
4.31 (iii) and the formula (|4.47p show that 

From and (|OT|) . we see that 

\cl,{z,z)\^\^{z,z/)\<\z% 

bmce z(x°,0) = 0, it follows that 

(j>{z{x'^,0),z{x°,0)) =0. (4.50) 

Using (|4.4ip again, we have that 



This together with (|4.26p and (|4.44p shows that 

4)iz{x°,0),z{x°,0)) = 32 
From (|4.50p and (|4.5ip . we see that 

^{x°,0) = I 32 



< 



a5(/)(i, Q)dT. 



(4.51) 



a50(r,O)dT,O 



□ 



Remark 4.7. From Lemma 14.61 it follows that the vacuum boundary Ttq in the coordinates (Sf,x^) is 
given by 

f / rS^/o 9i2(Cs",0)d?+zo ^ 

i'(rTj= <^ 32 / dS{T,Q)dT,{) 



< x° < Toj . 

We also note that the bounds on diz and 9^0(r, ^) from Proposition 14.31 (111) and (|4.42p . respectively, 
imply that the map 

fx'' So 9iz{(,x\o)di+zo 

32 / d^^{T, Q)dT 

Jo 

is strictly increasing. This shows that ^'(rj'Q) is just a reparameterization of Ttq, and that 

^To aiz(iTo,0)di+zo 



0<x^ < 



a^(/)(T,o)dT 



This is consistent with our assertion from Section r4.1l that the vacuum boundary in the (a;°, x^) coordinates 
is contained in the line x^ — 0. This fact combined with w = ^J^^do (see ()4.2p ) shows that the vacuum 
boundary moves with the fluid as noted previously. 
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4.5. Solutions with non-zero acceleration at the vacuum boundary. With the vahdity of the 
coordinate transformation ()4.12p established, we now turn to showing the that the maps {z, z,ip, cj)} 
determine a solution to the EFW equations that have non-zero acceleration at the boundary. We begin 
by letting 

g = ^/*g = g^^ dx^^ dx" and Hj = = 

denote the metric and the Frauendiener- Walton vector field. By (|4.2|) and (|4.7|) . the coordinate compo- 
nents of g and w are given by 

• (l + cos(z))^ 
{g,u) 4 I (4-52) 



(1 - cos(z))2 



and 



V2- 



where (see Lemma |4. 6 



- J^' [ ) (4.53) 

\d^(l){z,z) dg(l){z,z)J ^ ^ 

We also observe that the norm w^, the square of the sound speed s^, and the proper energy density p are 
easily computed to be 

= ii^^, (4.55) 
s'-2( '—^) , (4.56) 



1 + COs(2) 



and 



1 1 - cos(z) ,^ ^ , 

using the formulas (j^ . and (g^. From (11351) and we then obtain the following 

formula for the fluid two-velocity 

The analysis contained in Sections [3| 14.11 1421 14.31 and 14.41 guarantee that the pair {^^^i/, w"^} defined 
by the formulas ()4.52p . (|4.53p . and ()4.54p determine a solution to the EFW equations (|2.ip on the 
spacetime region Uto,s- This, in turn, shows that {g^,^, p}, with p and given by the formulas (|4.57p 
and (|4.58p . is a solution to the Euler equations (|1.2p - (ll.3p on Utq^s- 

With existence established, we are left with calculating the norm of the fiuid acceleration at the 
boundary. We begin by observing that 

(1 -|-cos(z))4 /I 

"32 (l-cos(z)) 

and 



9 = - ^TT^^dV'dV'-T^ r::7^d0d./) (4.59) 



_b y2(l + cos(z))4 



= dV' (4.60) 
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where = QfivW^ dx" . Using (|4.34p . we see that 



and hence by (I4.39p . that 



dijj = ^^_^ d^4>{z, z^) dz - 



sin'^(z) 



4z3 dt(t>{z, 7^ 
sm^(z) z"* 



■dz' 



Also, similar calculations show that 

-^dq 
z^ 

Next, we introduce the dual basis 



dt(j){z, 



■dz + 29t0(z,z'^)dz2. 



(4.61) 



(4.62) 



(4.63) 



f = dz and = dz^ 



(4.64) 

As we shall see, the components of the metric and Fraueniener- Walton covector field with respect to this 
frame have finite limits at the vacuum boundary even though some of the {x^) coordinates components 
are singular there. It is worthwhile noting that since this basis arises from the coordinate^ (z,z^), we 
could have introduced yet one more coordinate transformation to investigate the regularity of the fields 
{g, w} near the vacuum boundary. However, it is simpler just to work with the basis (|4.64p without 
introducing another explicit coordinate transformation. 
Writing the metric g and co-vector field tJ}' as 

g = (jijO^O^ and ^ Wi9\ 

a straightforward calculation using (|4.59p - (|4.64p shows that 

512(950(z(xO,O),O))^ 64a50(z(xO,O),O)a2^0(z(xO,O),O) 
64a5^(z(a;0, 0), 0)a2^^(z(a;0, 0), 0) -32(55^(z(xO, 0), 0))' + 8(92^^(z(x0, 0), 0))' 



(f^jlr^o) 
and 



(iD.lr^J = (329^0(2(^0, 0),0) 4a2^0(z(s",O),O) 



From (|4.65p . we find that 



(4.65) 
(4.66) 



- To J 



( 4(O;0(z(s°,O),O))'-(3,^g0(z(g°,O),O))' Of;0(z(5°,O),O) ^ 

2O48(a{0(z(2O,O),O)) 256(ac<^(z(sO,0),0))' 



afj<^(z(2<',o),o) 

256(ae<^(z(s0,0),0))^ 



32(9{0(z(£SO,O),O)) / 



Using this and (I4.66p . we get that 



V2 



(w^'Ir^J = I 3294^(z(iO,0),0) 
which we can use to write the frame components of h (see (|1.4p ) as 

I (o,";0(z(g°,O),O))' 3|g,^(z(a",0),0) ^ 

3f;0(z(g°,O),O) 1 

V 256(a{0(z(2O,o),O))' 32(9j0(z(2O,O),O)) / 



(4.67) 



^^Recall that it was shown in Section [4?4l that (z, z) satisfies dz Adz = {{diz)^ — (Sqz)-^) dx" Ada;-"^ with (((9iz)^ — (9o^)^) 
non-vanishing on Utq,&- This show that {z,z) define a coordinate system and it follows that does also. 
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Letting 

we get from (|4.56p that 
From this, we see that 



sin(z) 



z (l+cos(z))2 



(dsllr.J 



which, with the help of (|4.67p . shows that 



1 



Ids'l^lr^ = 3 T. 

" 128(a4</)(z(so,0),0)) 

This and p.lOp allows us to conclude that the norm of the fluid acceleration at the vacuum boundary is 
non-zero and is given by the formula 

|ak|r,„ = (4.68) 



/T289^0(z(iO,O),O) 
We summarize the above results in the following Theorem. 

Theorem 4.8. Suppose k = 3, z is the solution to the wave equation (j4.20l) from Provosition lA/dl z is a 
defined by (|4.26p , (f) and (j) are the maps from Proposition 14.51 and ip is given (I4.44p . Then for Tq and Sq 
small enough, the trivla^'^ {Stiu: PiV'^} determined by 



i9fii^) ~ J ill 

1 1 — cos(z) 

P = 7^ 

K 1 + cos(2) 



/[l + cos(z))4 

V 



32 




(1-C0S(Z))2, 



and 



where 



'J \Ti 



(l-cos(z))4 * V 

' d g(f){z, z) dz(l){z,z ) ' 
P{z) _ 

dz(f>{z,z) di(j){z,z) 



diz doz 
doz diz 



defines a solution of the Euler equations (jl.2p - (ll.3p on the spacetime region Uto,S- Moreover, the fluid 
acceleration |a|^ is non-zero on the vacuum boundary and given by the formula 

\o.\h\rT = . — 

" VT28d^(j){z{xO,0),0) 

/or < x° < Tq. 

^■^Here, the p, v'^, and are given in the (5;'^) coordinates with the Jacobian matrix arising from the transformation 
from the to the (x^^) coordinates. In the {x'^) coordinates, p, C^*, and are determined by 112.611 . 114.211 . and I I4.7II . 
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5. Exact Solutions 

In this section, we show that it is possible to determine certain classes of solutions to the Euler 
equations with non-zero acceleration at the boundary that are exact in the sense they are determined 
up to quadrature. First, we note that z and z given by the formulas (I4.25|) and (|4.26|) . respectively, are 
determined up to integrals. Next, by inspection, we observe that 

(/)(z, z) = (ci + C2z) (j^-^ cos(z) + ^ cos3(z)^ (5.1) 

satisfies the wave equation (|4.35p . and in fact, satisfies the conditions (|4.36p . (|4.37p . and (|4.38p . In this 
case, the map is given by 



cos(Cl/4) + J_ cos3(el/4)^ , 

Vl2 8 '^ 24 7' 



4>{t, = (ci + C2t)) - - - cos(Ci/4) + cos^e^^) , (5.2) 



and it is not difficult using the powerseries expansion for cos(x) that (j){t, ^) admits the follow expansion 
about ^ = 

= CO (^«-<^«"^H-^«»)+0«='»). (,3) 

Next, using ()4.43p . it follows from ()5.ip that the map ip is given by 

/ 1 1 cos(z) ln(csc(£) - cot(z)) 1 
■0(z, z) =C2 I -- + - ln(2) - -^^-^-^ 1 i J 



3 3 ' ' 3sin^(z) 3 2sin''(z) 

cot^(z) ln(sin(z))\ / C2 2 



, , ,ciz + — . (5.4) 
6 3 J \ 2 J ^ ' 

Together, {z{x),z{x)} determined by (|4.25p and (|4.26l) . and {<j),4>,ip} determined by ((5?T|) . (fO)) . and 
(|5.4p specify completely a solution to the Euler equations via the formulas in Theorem 14.81 Moreover, we 
see from ()4.68p and (|5.3p that the acceleration at the vaccuum boundary for these solutions is given by 
the formula 

l-l I 32 

yT28(ci+C2z(xO,0))' 

and in particular, is constant if C2 = and time varying otherwise. 
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Appendix A. The singular wave equation + = 

Our goal in this Appendix is to prove the existence and regularity of solutions to the singular, linear 
wave equation (I4.35p . To simplify notation, we will set (z, z) = (t, x) as we will be thinking here of z and 
z as a time and space coordinate, respectively. Letting 



^/^^{x)' 

we see from (|4.35p that $ satisfies the equation 

d^^ + H^ = (A.l) 
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where H is the operator 

1 d 



A.l. The Friedrichs extension of H. First, we observe that a simple integration by parts argument 
shows that H defined (jA.2|) is symmetric on the domain C^(0,7r/2) with respect to the standard 
inner product 

Jo 

We also observe, again by integration by parts, that 

H = L^L 



where 



1 



and f is the adjoint given explicitly by 



1 



This shows that H is a non-negative symmetric operator and the quadratic form 



'?if(4'l,$2) = (i$l|L$2) 

= (A.4) 
||<i>||?,. -||i<i>||i. + ||<i>lli., (A.5) 



associated to H satisfies 
and in particular. 
Defining the norm 
we let 

denote the completion of C^{0,tt/2) with respect to the norm (|A.5p . Then by Theorem X.23 of [17] . 
the self-adjoint Friedrichs extension of H, which we also denote by H, exists and is defined on a dense 
domain 

D{H) CH^ C L2(0,7r/2). 

A. 2. Existence. In order to discuss existence for the wave equation, we need define the norms 

fc 

= Il^'^'^lli^ (-^'^' - (H^^y ) (A.6) 

3=0 

and the spaces 

k 

n'^^f] D{W/^) C L2(0, 7r/2), (A.7) 
j=o 

where H^/"^ is the positive square root of H. 

The following Theorem which guarantees existence and uniqueness follows from well-known results 
existence and uniqueness results for abstract wave equations. See [HI HH] for details. 
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Theorem A.l. Suppose fc e No and (<i>o,$i) G "H + x Ti . Then there exists a unique solution 
$ e n*^±oC^([0,oo),H''+^"^) to t/ie mifia; value problem 

$|t=o, -^l*=o 1 = ($0,$l)- 

A. 3. Regularity near the boundary. Away from a; = 0, the norms || • \\-j^k are equivalent to the 
standard Sobolev norms and, hence, any solution to (|A.ip in "H*^ will also lie in Hf^^^{0,TT/2). However, 
the norms || • Wy^k are not uniformly equivalent to the standard Sobolev norms as x approaches zero. 
Consequently, some work is needed to determined the regularity of functions lying in these spaces for 
small X. Our main tools to establish the regularity will be the following: 

(i) Sobolev's inequality: Suppose sp < 1, ,s G M, and p G (1, oo). Then 

ll*llLP/(i-=p)(a,6) < ll$lk=.P(a,;,)- (A.8) 

for all $ G H'''P{a, b). Here, H'''P{a, b) denotes the fractional Sobolev spaces which coincide with 
standard ones for s G N. We employ the usual notation H^{a^ b) = _ff^'^(a, b), and we note that 
for < s < 1 the fractional norm ||$||_f/s(a,6) can be written as 

mhia.b) = mhia.b) + f f^-^^^—§^dxdy. (A.9) 

Ja Ja P y\ 

(iii) Morrey's inequality Suppose p G (l,oo]. Then 

ll*llco.i-i/p < ||$||ifi-(a,;,) (A.IO) 

for aU $ G H^'P{a,b). 

(ii) Fractional order weighted Hardy's inequality: For < d < 1, the following inequality 



\^\\H^-''{Q,b) ^ II*IIl2(o,6) 



x'' — 

dx 



L2(o,6) 



(A.ll) 



follows directly from the fractional order weighted Hardy's inequality [131 Theorem 5.3] and the 
definition of the fractional norm (|A.9p . We note that this type of Hardy inequality was also used 
in the existence results of [H El jTll [12] . 

To begin, we introduce a new coordinate 

^ = [ fi{x)dx, (A.12) 
Jo 

and let 

Co = / M(s)ds. 

"'0 

We use the notation L| to denote the space with respect to the coordinate ^ on the interval (0,^0)7 
or equivalently, the measure 

d^ = fi{x)dx (A.13) 
on the interval (0,7r/2). We use the following notation for the norm: 

ii'i>iiL^ = ^^"<i'm- 
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We will also use the notation i/| when referring to the Sobolev spaces with respect to the variable ^, 
or equivalently, with the spaces defined using the differential operator 

d \ d , . s 

(A.f4) 



fi{x) dx 

and the measure jXTSj. In the following, we use the same notation to denote a function whether thought 
of as a function of a; or a function of ^ whenever this does not lead to an ambiguity. 
Since is analytic, > 0, and ^{x) — 0{x'^), it follows that 

^ - a; 



(A.15) 



for X near zero, or equivalently 

for ^ near zero. Thus, in particular, 

for ^ near zero. 

Lemma A. 2. Suppose ^ e H'^ . Then 



X ~ 



(A.16) 



Proof. First, we observe that the identities 



m^dx = 



\LWdx = 



|(V7i<J>) 



d^, 



and 



d^ 



(A.17) 
(A.18) 

(A.19) 



follo w dire ctly from the definitions (|A.2p . (|A.3[) . (|A.12p . (|A.13p . and (|A.14p . Integrating (|A.17p . (|A.18[) . 
and (|A.19P then gives 



<ii'i>ir+iiL<i>iii. 



(A.20) 



Since 



by (|A.4p . the proof follows from (|A.20p . 
Lemma A. 3. Suppose $ e H''. Then 

and 



mil. ^ \\H'/H\\i. 



□ 
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Proof. First, we observe that 



^3/4 d^2 

Next, we see from (|A.2p that 



< 



< 



^1 



(by (KM) 

(by Lemma IA.2P 



(A.21) 



/ ^d^ 



de 



de 



and so upon integrating, we find, after using (|A.16p . that 



de 



< 



2 



Setting 



the two inequahties (|A.2ip and (|A.22p show that 



d^ 



^3/4 



^3/4 



de 



< ml.. 



Letting 
we compute 



sup 1/(01 

A/2<C<A 



sup I /a (01 

1/2<C<1 



< 



< 



< 



ll/A||ffi(l, 



2) 



mei'de 



1/2 



/" 


dVx 


Jl/2 


de 



d^ 



1/2 



^3/2 



1/2 



de 



(AO 



d^ 



J\/2 



d'f 



2(0 



d^ 



(A.22) 



(A.23) 



(by (KM) 



jx/2 ' 

<Ai/2||$||?,., 

where in deriving the lass inequahty we used (jA.23p ). From this, we conclude that 

1/(01 <ll<i>ll«^IO'/'- 

In particular, this shows that 

lim /(O = 0. 

Next, we see that from (|A.23P and the inequality (jA.lip that 

ll/IUv4 < ||a>ll«4. 

But since 

l|/l!c°^3/4 < 11/11^1,4 < ||/|1^V4 



(A.24) 
(A.25) 
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by (jA.sp and ()A.10| . we see that 



25 



This together with (|A.25[) shows that 
which in turn imphes that 



< 



< 



2-e 



for any e > 0. Integrating, we find that 



(A.26) 



and this inequahty combined with Lemma IA.2I and the inequahty (|A.11[) shows that 



Applying Morrey's inequahty (jA.lOp . we arrive at 

d 



< ll$l 



-^0,1/4 



Next, we observe that 



+ 



i 5 

by Lemma IA.2I and (jA.I6| . The same argument used to derive (|A.24p from (jA.23| shows that 

\Vmm\<\mn-^e^'- (A.27) 

Together, Lemma [A. 2[ the inequahty (jA.26|) . and (jA.lip show that 
and hence that 

||V^a>|| < ||$||„4 (A.28) 

by Morrey's inequahty (jA.10[) . Finally, since 



hm vmm = 



by (|X27)) . we get from (|X28|) that 



□ 



Theorem A. 4. Suppose ($o,<&i) eV.^ xH^ and $ e Hj^QC^ {[O^oo),^^ ^) is the solution to the wave 
equation (jA.ip from Theorem \A.l\ Then there exists a map (f){t,£,) that satisfies the following: 
(i) 

4> e C\[0, oo), Ci^i/2(o, (^/2)4)) n f] ([0, oo), C^-^{0, (7r/2)4)), 

]=0 
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(ii) for any t £ (0, oo), 

\dtkt.o\ + m.o\<^ 

forO<C< (7r/2)4, and 

(iii) 



for all {t,x) e (0,oo) x (0,7r/2). 
Moreover, if there exists a positive constant c such that 

^54VMaO*(0,2;)) > c> 

for all X G (0,7r/2), t/ien f/iere exists aT > such that 

/or all {t,x) e (0,r) X (0,(7r/2)4). 

Proof. Statements (i)-(iii) follows directly from the regularity statement $ e n^^QC-'([0,(X)),'H^~^), 
Lemma (jA.Sp . and that fact that C H^^^ which follows from Morrey's inequality and the inclu- 
sion C H^^^{0, tt). Note that we are also using the fact that the variable ^ is uniformly equivalent to 
x"^, see (|A.15I) . 



For the final statement, suppose that 

1 



for X e (0,7r/2). Then differentiating ^ fi{x)^(t, x) = (j){t,x'^) gives 

4a:" 



d^{^Ji{x)'^{Q,x))>c>0 (A.29) 



/i(a;) ' 
Since there exists a non-zero constant C such that 

1 

for < x < 7r/2, it follows from (|A.29p and (|A.30p that 



5x(v^*(i,a:)) = —-d^ct>{t,x^). (A.30) 



c 



for all X G (0, 7r/2), or equivalently 



5^0(0, > ^ >0 



5e0(O,C)>^>O 



for all ^ e (0, (7r/2)^). From the continuity of d^(j){t,^), we see that there exists a T > such that 

for all (t, G [0, T) x (0, (7r/2)4). □ 



SOLUTIONS TO THE RELATIVISTIC EULER EQUATIONS WITH A VACUUM BOUNDARY 



27 



References 

[1] L.C. Evans, Partial Differential Equations, AMS, 1998 

[2] Y. Choquet-Bruhat, C. De Witt-Morette, and M. Dillard-Bleick, Analysis, Manifolds, and Physics Part I: Basics, 
revised ed., North-Holland, Amsterdam, 1996 

[3] D. Christodoulou, The Formation of Shocks in 3- Dimensional Fluids, EMS, 2007 

[4] D. Coutand, H. Lindblad, and S. SlikoUcr, A priori estimates for the free-boundary 3-D compressible Euler equations 

in physical vacuum, Comm. Math. Phys. 296 (2010), 
[5] D. Coutand and S. Shkoller, Well-posedness in smooth function spaces for the moving-boundary 1-D compressible 

Euler equations in physical vacuum, preprint [aj:Xiv:0910.3136] 
[6] D. Coutand and S. Shkoller, Well-posedness in smooth function spaces for the moving-boundary 3-d compressible Euler 

equations in physical vacuum, preprint [arXiv:1003.4721] 
[7] J. Frauendiener, A note on the relativistic Euler equations. Class. Quantum Grav. 20 (2003), L193-L196 
[8] H. Friedrich, Evolution equations for gravitating ideal fluid bodies in general relativity Phys. Rev. D 57 (1998), 2317- 

2322 

[9] J. A. Goldstein, Semigroups of Linear Operators and Applications, Oxford University Press, 1985 
[10] J. A. Goldstein adn M. Wacker, The energy space and norm growth for abstract wave equations, App. Math. Lett. 16 
(2003), 767-772 

[11] J. Jang and N. Masmoudi, Well-posedness for compressible Euler with physical vacuum singularity. Comm. Pure Appl. 
Math. 62 (2009), 1327-1385 

[12] J. Jang and N. Masmoudi, Well-posedness of compressible Euler equations in a physical vacuum, preprint 
[arXiv: 1005.4441] 

[13] A. Kufner and L.E. Persson, Weighted Ineqaulities of Hardy Type, World Scientific, 2003 
[14] S. Lang, Fundamentals of Differential Geometry, Springer, 1999 

[15] P.G. LcFloch and S. Ukai, A symmetrization of the relativistic Euler equations in several spatial variables. Kinetic 

and Related Models 2 (2009), 275-292 
[16] T. Makino, "On a local existence theorem for the evolution equation of gaseous stars" , in Patterns and Waves, edited 

by T. Nishida, M. Mimura, and H. Fujii, North-Holland, Amsterdam, 1986 
[17] M. Reed and B. Simon, Methods of modern mathematical physics II: Fourier analysis, self-adjointness. Academic 

Press, 1975 

[18] A.D. Rendall, The initial value problem for a class of general relativistic fluid bodes, J. Math. Phys. 33 (1992), 
1047-1053 

[19] M. Visscr and C. Molina-Pan's, Acoustic geometry for general relativistic barotropic irrotational fluid flow. New Journal 

of Physics 12 (2010), 095014 (18pp) 
[20] R.M. Wald, GeneraZ Relativity, University of Chicago Press, 1984 

[21] R.A. Walton, A symmetric hyperbolic structure for isentropic relativistic perfect fluids, Houston J. Math. 31 (2005), 
145-160 

School of Mathematical Sciences, Monash University, VIC 3800, Australia 
E-mail address: todd.oliynykesci.monash.edu.au 



